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ABSTRACT 

Aims. We observed the high velocity outflow originating from NGC6334 I in several CO transitions with the APEX telescope, with the goal of 
deriving the physical parameters of the gas. 

Methods. Using an LVG analysis, we studied line ratios between the CO(3-2), CO(4-3), and CO(7-6) data as a function of the density and of 
the kinetic temperature of the gas. An upper limit on the CO column density is derived by comparison with 13 CO data. 

Results. We constrained the temperature to be higher than 50 K and the IF density to values higher than n ~ 10 4 cirr 3 towards the peak 
position in the red lobe, while T > 15 K and n > 10 3 cirr 3 are derived towards the peak position in the blue lobe. The physical parameters of 
the outflow, its mass and its energetics, have been computed using the temperatures derived from this analysis. 

Conclusions. We conclude that high kinetic temperatures are present in the outflow and traced by high excitation CO lines. Observations of 
high-J CO lines are thus needed to infer reliable values of the kinetic temperatures and of the other physical parameters in outflows. 
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1. Introduction 

Despite molecular outflows being a well-studied phenomenon 
associated with young stellar objects of all masses, many ques- 
tions on their driving mechanisms and on their interaction with 
the surrounding material are still unanswered. The case of mas- 
sive outflows, which only recently have been investigated at 
high spatial resolution (e.g. lBeuther et alJl2002at Isollins et alJ 
2004), is even more problematic. Usually, low rotational (J) 
CO transitions are used to trace molecular outflows and derive 
their physical parameters and energetics. Howe ver, these lines 
trace o nly cold gas, while high-J transitions ("e.g. lBeuther et alJ 
2002a) have also been detected, clearly pointing to hot, dense 
gas associated with the outflow. High resolution observations 
of high-J CO transitions are needed to get a better understand- 
ing of the phenomenon. 

The region NGC6334, located at a distance of 1.7 kpc from 
the Sun (Neckel 1978), harbours a series of luminous star- 
forming regions at various stages of evolution. NGC6 334 I is 
the d ominant source at (sub)millimeter wavelengths (Sandel 
2000), with a compact Hn region at its c enter ( Carral et al 
19971) Sub-mm line surveys JSandelil994llMcCutcheon et al 



20001 iThorwirth et alJEo03l Schilke et al. this volume) reveal 
rich molecular line emission from a dens e, hot core of size 
< 10". Two outflows originate from it JBachiller & Cernic haro 
ll990tlMcCutcheon et all2000l) . 

In this Letter, we present high frequency 13 CO and CO ob- 
servations of the main outflow originating from NGC6334 I. 



The physical parameters of the outflow are derived by a Large 
Velocity Gradient (LVG) analysis of several CO transitions. 

2. Observations 

The observations of NGC6334 I (a J2 ooo =17:20:53.35, 
£j2000 =-35:47:01.50) were performed with the Atacama 
Pathfinder Experiment 12 m telescope (APEX 1 , Giisten et al. 
2006, this volume) during June, September, and November, 
2005. The APEX-2A (Risacher et al. 2006, this volume) and 
FLASH (Heyminck et al. 2006, this volume) receivers were 
used to map the source in CO(3-2) and CO(4-3). The sec- 
ond channel of the FLASH receiver was tuned to CO(7-6); 
we performed long integration observations of the two posi- 
tions where the blue and redshifted emission in CO(4-3) is 
strongest, at the offset positions (-10", -14") and (17.5", 14"). 
Additional observations of the same two positions in 13 CO(3- 
2) were performed in November 2005 with the APEX-2A re- 
ceiver. All observations were taken in the lower side band. For 
the CO(4-3) and CO(7-6) observations, the system tempera- 
tures are around 800 K and 3400 K, at 461 GHz and 807 GHz 
respectively, while the system temperatures of the APEX-2A 
receiver are of the order of 260 K for CO(3-2) and 480 K 
for 13 CO(3-2). The system temperatures of the 13 CO(3-2) are 
higher, as they were performed in November 2005 when the 
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source was already in the Sun avoidance limit and was ob- 
served at low declination (< 20°) after sunset. Pointing was 
checked on the continuum of NGC6334 I itself, during the 
CO(4-3) observations, and was found to be accurate to ~ 5". 
On the other hand, problems with the pointing were encoun- 
tered during the CO(3-2) run and the data show a system- 
atic shift of ~ 8" with respect to the CO(4-3) observations. 
Alignment of the CO(3-2) map to the CO(4-3) is possible since 
several molecules are detected in the two bands that have high 
excitation conditions and peak on the hot core. We estimate the 
accuracy to be good to 2-3", since it was based on the align- 
ment of hot core lines, which trace the (within our beam) point 
like hot core. 

The FFTS spectrometer (Klein et al. 2006, this volume) was 
used with resolutions of 0.1 km s" 1 for CO(3-2), 0.3 km s" 1 for 
CO(4-3) and CO(7-6), and 0.4 km s" 1 for 13 CO(3-2). The cali- 
bration was performed by using the APECS software (Muders 
et al. 2006, this volume). Antenna temperatures were con- 
verted to main-beam temperatures by using forward and beam 
efficiencies of 0.97 and 0.74 at 345 GHz, 0.95 and 0.60 at 
464 GHz, and 0.95 and 0.43 at 809 GHz. 

The CO(4-3) map is shown in Fig.^ while spectra towards 
the (-10 ", -14") and (17.5", 14") posit ions are presented in 
Figs. 12131 Bachill er & Cernicharol Jl990l) detected emission at 
-100 km s" 1 in the CO(2-l) spectra; they did not identify it 
with a CO velocity component, even though they could not 
assign the feature to any other molecule, because it peaks at 
the hot core position. The detection of a similar feature at - 
100 km s" 1 in our CO(3-2) spectra (Fig. [3ji supports the in- 
terpretation that this is indeed another CO high velocity com- 
ponent. Alternatively, the feature can be identified as CH3OH 
(I61 — > 152, Ei ow ~ 316 K) as it peaks on the hot core. A 
possible interpretation is that this feature is a CO component, 
since high excitation methanol lines are unlikely to be detected 
at positions distant from hot cores, and that contamination with 
C H3OH happens at the hot co re position. 

iMcCutcheon et alJ ( 120001) show that the absorption at 
6.5 km s" 1 (see Figs. |3-|3|l is real, not coming from the off 
position used for position switching. 

3. Physical parameters of the outflow 

To determine the physical parameters in the outflow, we used 
the CO observations and a spherically symmetric LVG statisti- 
cal equilibrium code, with the cosmic background as the only 
radiation field. To compare the observations at the three differ- 
ent frequencies with the LVG results, we smoothed the CO(4-3) 
and CO(7-6) data to the same spatial resolution of the CO(3-2), 
which is 18". To reduce the noise level in the spectra, we re- 
sampled them to a resolution of 5 km s" 1 . 
The non-detection of the 13 CO(3-2) transition in the wings 
gives an upper limit to the CO column density of 10 17 cm" 2 
for the (-10", -14") and (17.5", 14") positions. In this cal- 
culation, we used the relation between 12 C and 13 C given by 
Wilson & Roo dllll994h and a galactocentric distance of 6.8 kpc 
( Kra emer et alJl998l) . All models were run for this value of the 
column densiy. At this column density, the CO transitions are 
optically thin, with the exception of CO(3-2), which is opti- 
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Fig. 1. Map of the CO(4-3) emission in the 7 < v < 55 km/s 
(red, dashed line) and -78 < v < -15 km/s (blue, solid line) 
ranges. Contours are from 100 to 600 K km/s by steps of 100. 
The star marks the NGC6334 I position; the beam (~ 13.5") is 
indicated in the top right corner. 
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Fig. 2. CO(J=3,4,7) and 13 CO(3-2) spectra at the offset 
(17.5", 14") from the reference position (anooo =17:20:53.35, 
£12000 = -35 : 47 : 01.50). The CO(7-6) and CO(4-3) spectra 
are smoothed to the CO(3-2) spatial resolution. 



cally thick for T< 50 K, n(H2) < 10 4 cm" 3 . Hence, our results 
should not depend much on the column density used in the cal- 
culations. The code inputs are the density and the kinetic tem- 
perature of the gas. Since a direct estimate of the latter is not 
available, we ran models for several densities and temperatures. 

Source sizes play an important role in comparing the obser- 
vations to the results of the LVG simulations, as the main beam 
temperature of a transition is equal to sl s ~ Qi x T/, where s and Q 
are the source size and the beam size in arcsec, and T/ is the line 
intensity from the LVG analysis. In our case, the beam is the 
same for all the transitions, but the source sizes may differ as 
the three lines have different excitation conditions. Therefore, 
line ratios may also be affected by the uncertainties in the beam 
filling factors. In the following, we will analyse the line ratios 
between the three transitions (upper J divided by lower J), as- 
suming that all transitions are emitted by the same volume of 
gas. These are lower limits to the effective line ratios, as less 
excited transitions are likely to come from more extended re- 
gions. 
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Fig. 3. CO(J=3,4,7) and 13 CO(3-2) spectra at the 
offset (-10", -14") from the reference position 
(aj2000 =17:20:53.35, 5 J20 oo = -35 : 47 : 01.50). The 
CO(7-6) and CO(4-3) spectra are smoothed to the CO(3-2) 
spatial resolution. 



In Fig. E] the three different line ratios at (17.5", 14"), 
Vhr - 32 km s~ , are plotted as a function of density and ki- 
netic temperature. Since we have ignored the effects of beam 
dilution, which are likely to be important at least for the CO(7- 
6)/CO(3-2) line ratio, we can only estimate lower limits to the 
gas parameters. We find that the temperature must be higher 
than 50 K, while the density is of the order of 10 4 crrT 3 , equiv- 
alently a thermal pressure of 5 x 10 5 crrT 3 K. Moreover, as the 
line ratios between the CO(7-6) line and the other transitions 
show, there is also a denser, warmer component, from which 
the CO(7-6) line is emitted. As already discussed, no correc- 
tions for the different beam filling factors are applied, and the 
values we used for the line ratios are likely to be lower lim- 
its. Therefore, the values of n — 7 x 10 4 cirT 3 and T — 70 K 
(corresponding to P thermal ~ 5 x 10 6 cm -3 K) that we derived 
for the component traced by CO(7-6) are lower limits (con- 
sistent with other stu dies of high-J transitions in outflows, e.g. 
iBeuther et alJl2002al) . Taking into account a calibration error 
of 0.2 for each ratio, we find 30 < T < 55 K, 7 x 10 3 < 
n < 2 x 10 4 cm 4 . For the warmer component traced by CO(7- 
6)/CO(4-3) 55 < T < 95 K, 3 x 10 4 < n < 1 x 10 5 cnr 3 . The 
results for the other velocity channels are very similar, as the 
line ratios between the different CO transitions do not signifi- 
cantly change at this position (F"ig. |2j. 

For the analysis on the (-10", -14") position, we only used 
the CO(4-3)/CO(3-2) ratio, since the signal-to-noise ratio in the 
CO(7-6) spectrum is not very high in the wings (~ 3). The 
lower limits on the temperature and the density of the gas vary 
along the wings, ranging between 15 and 50 K, 2 x 10 3 and 
2 X 10 4 crrT 3 (see Figs . |5 (a jf |5 (b) | > . Taking into account the un- 
certainties on the calibration, the lower limits we derive range 
between 10 and 30 K and n < 10 4 cirT 3 for Vi sr = -48 km s , 
and 20 < T < 56 K and 4 x 10 3 < n < 3 x 10 4 cm 3 for 
= -5810ns" 1 . 



Vis, 

The difference in temperatures and densities derived for the 
red and blue lobes may be caused by the (pre-existing) den- 
sity structure of the molecular cloud into which the outflow 
impinges and with which it interacts. Continuum maps of the 
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Fig. 4. Line ratios as a function of temperature and density to- 
ward the (17.5", 14") position. 
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Fig. 5. CO(4-3)/CO(3-2) line ratio as a function of temperature 
and density for two velocity channels toward the (-10", -14") 
position. 



region ISanded koOO) support this scenario, since the contin- 
uum emission is extended to the north of NGC6334 1. 

The CO column density of the high velocity gas can be de- 
rived from the integrated intensity, by assuming that the level 
populations are in LTE at the kinetic temperature derived in 
our calculations (30 K for the blue lobe, 50 K for the red one). 
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To avoid contamination from the core emission, we integrate 
between 10 and 59 km s for the redshifted gas, and -88 to 
-22 km s _1 for the blueshifted gas, which results in lower lim- 
its to the derived masses. The momentum p, the energy E, and 
the characteristic time scale t can be derived assuming that the 
outflow is a constant velocity flow wi th true velocity equal to 
the maximum velocity observed (e.g. Lada 1985; Beuther et al. 
2002b). Finally, the mass entrainment rate M to t, the mechani- 
cal force F, and the mechanical luminosity L can be derived 
from the other parameters. Assuming that the inclination angle 
between the outflow and th e line of s ight is 9, the following 
relations are used (from Beu ther et all2002bl but corrected for 
the inclination angle 6): 



N(H 2 ) = 10 



4 6.97 x 10 15 {F/T C 

— 2 7^™ j T co( 4-3)t 



M tot — (Nb X area/, + N r X area,) mu 

cosQ 



P = {M b X Vmax b + M r X V maXr ) 

cost) 



1/1 ,1 7 

E = ~, I ^M b x v L maXh + -M r x vf r 

cos 1 - 



t = 



cosQ 



IsinQ 



+ v„. 



(1) 

(2) 
(3) 
(4) 
(5) 



the analysis of CO(2-l) and C O(3-2) observations . Both 
IBachiller & Cernicharol Jl990l) and iMcCutcheon etalJ §000) 
do not apply any corrections for the inclination angle to their 
values. In Tabled the parameters are given for no corrections 
for the inclination angle and for 8 = 45°and 70°, as the small 
overlap between the redshifted and the blue emission and the 
emission at high velocity indicate tha t 9 is closer to 90° than 
to 0°. For comparison, results from IBachiller & Cernic haro 
(1990) and Mc Cutcheon et alJ d2000l) are also given. Our val- 
ues are in good agreement with previous studies; however, 
the mass we derive, and therefore other paramet e rs also , are 
smaller than what both IBachiller & Cernicharol dl990l) and 
McCutcheon et al. (2000) had derived, due to the higher ki- 
netic temperature we use for the red lobe. Morever, high- J tran- 
sitions trace only the warm gas, while lo w-J lines also trace 
a cold er component. Both this study and IMcCutcheon et all 
(2000) make use of the maximum observed velocity in 
the outflow to derive the en ergetics of the outflow, while 
IBachiller & Cernicharol Jl990h use the mean velocity. The es- 
timate of the mass is also affected by other uncertainties; the 
observed area of the outflow depends on the transitions used in 
the calculations and the excitation temperat ure used in Eq. (EJ 
has been derived for one position per lobe. Gibb et aD (|2003) 
also found lower outflow masses in CO(3-2) than in CO(1-0); 
ICabrit & Bertoutl (ri990) have shown that masses derived using 
Eq. (|2j are upper limits to their true value. 
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Table 1. Physical parameters of the outflow 







1.9 x 10 21 cm 


-2 






1.5 x 10 21 cm" 


-2 




no corr. 6 = 45° 


e = 70° 


M b [M ] 


1.1 


1.5 


3.2 


M r [M ] 


0.8 


1.2 


2.4 


M to , [M f 


1.9 


2.7 


5.6 


t [10 3 yr] 


4.4 


4.5 


1.6 


p [M Q km s' l ] b 


140 


200 


410 


£ [10 47 erg] 


1.0 


2.1 


8.9 


F m [M kms-' yr 1 ] 


0.03 


0.04 


0.3 


L m I /- ] 


780 


1550 


1800 


M fof [10- 3 M o yr-'] 


1.9 


2.6 


15 



(6) 4. Conclusions 
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Estimates of the 
flow from NGC6334 



physical parameters in the out- 
I have already been derived from 



We report observations in high excitation CO lines towards the 
high velocity outflow from NGC6334 1. By using an LVG anal- 
ysis, we have constrained the temperatures to be higher than 
50 K and the H2 density to values higher than n ~ 10 4 cm" 3 
towards the peak position in the red lobe, T > 15 K and 
n > 10 3 cirT 3 towards the peak position in the blue lobe. The 
physical parameters of the outflow, its mass and its energet- 
ics, have been computed using the temperature s derived from 
this a nalysis. Often (e.g. Beut her et al J l2002ht iKim & KurtH 
2006), low temperatures are assumed in computing outflow 
masses; our analysis has shown that high kinetic temperatures 
are present, and traced by high excitation CO lines, while they 
are not found when only low energy transitions are available. 
High spatial resolution maps of CO transitions at high-J are 
needed to infer reliable values of the kinetic temperatures and 
of the other physical parameters in the outflow. 
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